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Faraday Rotation and One-Way Propagation of Plasmon Waves 
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Abstract -We calculate the dispersion relations of plasmonic waves propagating along a chain 
of metallic nanoparticles in the presence of both a static magnetic field B and a liquid crystalline 
host. The dispersion relations are obtained using the quasistatic approximation and a dipole-dipole 
approximation to treat the interaction between surface plasmons on different nanoparticles. For 
a plasmons propagating along a particle chain in a nematic liquid crystalline host and a field 
parallel to the chain, we find a small, but finite, Faraday rotation angle. In a cholesteric liquid 
crystal host and an applied magnetic field parallel to the chain, the dispersion relations for left- 
and right-moving waves are found to be different. For some frequencies, the plasmonic wave may 
propagate only in one of the two directions. 
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Introduction. — Ordered arrays of metal particles 
jin dielectric hosts have many remarkable properties [HE]. 
For example, they support propagating modes which are 
linear superpositions of so-called ’’surface” or ’’particle” 
plasmons. In dilute suspensions of metallic nanoparticles, 
these surface plasmons give rise to characteristic absorp¬ 
tion peaks, in the near infrared or visible, which play an 
important role in their optical response [THU]. For ordered 
arrays, if both the particle dimensions and the interpar¬ 
ticle separation are much smaller than the wavelength of 
light, one can readily calculate the dispersion relations for 
both transverse (T) and longitudinal (L) waves propagat¬ 
ing along the chain, using the quasistatic approximation, 
in which the curl of the electric field is neglected. 


^ In a previous paper, we calculated these dispersion rela- 
5 ^ tions for metallic chains immersed in an anisotropic host, 
^such as a nematic or cholesteric liquid crystal (NLC or 
CLC) [S]. Here, we consider the additional effects of a 
static magnetic field applied either parallel and perpen¬ 
dicular to the chain. For a parallel magnetic field orienta¬ 
tion we find that a linearly polarized T wave undergoes a 
Faraday rotation as it propagates along the chain. For a 
field of 2 Tesla, and reasonable parameters for the metal, 
this Faraday rotation is around 1/2 a millidegree per ten 
interparticle spacings. These results suggest that either a 
nematic liquid crystal (NLC) host or an applied magnetic 
field could be used as an additional “control knob” to ma¬ 


nipulate the properties of the propagating waves in some 
desired way. 

We also consider the same problem when the host is a 
cholesteric liquid crystal (CLC). In this case, if the mag¬ 
netic field is parallel to the chain and the director rotates 
about the chain axis with a finite pitch angle, we show 
that the frequencies of left- and right-propagating waves 
are, in general, not equal. This difference opens up the 
possibility that, for certain frequencies, a linearly polar¬ 
ized wave can propagate along the chain only in a single 
direction. Thus, this geometry provides a possible real¬ 
ization of a one-way waveguide. This realization is quite 
different from other proposals for one-wave waveguiding, 
such as electromagnetic waves propagating at the inter¬ 
face between two different materials [231E3], along chains 
of ellipsoidal particles in a spiral configuration , and 

in arrays of gyromagnetic crystals m- 

The remainder of this article is organized as follows: 
First, we use the formalism of Ref. [5] to determine the 
dispersion relations for the L and T waves in the presence 
of an anisotropic host and a static magnetic field. Next, 
we give simple numerical examples and finally we provide 
a brief concluding discussion. 

Formalism. — We consider a chain of identical metal 
nanoparticles, each a sphere of radius a, arranged in a one¬ 
dimensional periodic lattice along the 2 axis, with lattice 
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spacing d, so that the particle is assumed centered at 
(0,0, nd) (—00 <n < + 00 ). The propagation of plasmonic 
waves along such a chain of nanoparticles has already been 
considered extensively for the case of isotropic metal par¬ 
ticles embedded in a homogeneous, isotropic medium m- 
In the present work, we calculate, within the quasistatic 
approximation, how the plasmon dispersion relations are 
modified when the particle chain is immersed in both an 
anisotropic dielectric, such as an NLC or CLC, and a static 
magnetic field. We thus generalize earlier work in which 
an anisotropic host is considered without the magnetic 
held 012]. 

In the absence of a magnetic held, the medium inside 
the metallic particles is assumed to have a scalar dielec¬ 
tric function. If there is a magnetic held B parallel to 
the chain (which we take to lie along the z axis), the di¬ 
electric function of the particles becomes a tensor, In 
the Drude approximation |28j . the diagonal components 
are e(a;), while Cxy = —Syx = iA(uj) and all other compo¬ 
nents vanish. In this case, the components of the dielectric 
tensor take the form 


and 


e(w) = 1 - 


CJ" 


' + i/r) 


1 - 


lj: 


A{uj) 


WpT uJcT 

W (1 — IUJtY 


( 1 ) 

( 2 ) 


where ujp is the plasma frequency, r is a relaxation time, 
and uJc is the cyclotron frequency, and the second limit 
is applicable when wt 00 . We will use Gaussian units 
throughout. 

The dielectric function of the liquid crystal host, for ei¬ 
ther the NLC or CLC case, is taken to be that described 
in Ref. [Sj The dispersion relations for the surface plas¬ 
mon waves are determined within the formalism of Ref. [Sj 
Specifically, we write down a set of self-consistent equa¬ 
tions for the coupled dipole moments; these are given in 
Ref. [5|as Eq. (9), and repeated here for reference: 
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Here ^ 

i = Se(i- f (4) 

is a “t-matrix” describing the scattering properties of the 
metallic sphere in the surrounding material, Q and F are 
a 3 X 3 Green’s function and depolarization matrix given 
in Ref. [5] 1 is the 3x3 identity matrix, and Se = — £h, 

where ih is the dielectric tensor of the liquid crystal host. 

Nematic Liquid Crystal. We first consider a chain of 
metallic particles placed in an NLC host with B||z and 
parallel to the liquid crystal director n. Using the for¬ 
malism of Ref. [5] combined with Eq. ©, we obtain two 
coupled sets of linear equations for the transverse (T) com¬ 
ponents of the polarization, pnx and Pny The solutions are 
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Fig. 1: Blue symbols (x’s and -l-’s): Dispersion relations for 
left and right circularly polarized T plasmon waves propagating 
along a chain of metallic nanoparticles immersed in a NLC with 
both the director and a magnetic field parallel to a chain. The 
blue symbols correspond to metal particles are described by a 
Drude dielectric function with LOpT = 100 and lUcj^p = 3.5 x 
10~®. The red symbols (open squares and triangles) represent 
the same dispersion relations as the blue symbols, but assuming 
no single-particle damping, corresponding to cupT = 00 . In both 
cases, the splitting between left and right circularly polarized 
waves is not visible on the scale of the figure. For Wp = 1.0 x 
10^® sec“^, the chosen ojc/wp corresponds to R ~ 2 Tesla. 
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found to be left- and right-circularly polarized transverse 
waves with frequency lu and wave number k±^ where the 
frequencies and wave numbers are connected by the dis¬ 
persion relations in the nearest-neighbor approximation 
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where we use the notation of Ref. [Sj The longitudinal 
{L or z) mode is unaffected by the magnetic held. Since 
the frequency-dependences of both e{uj) and A{uj) are as¬ 
sumed known, these equations represent implicit relations 
between w and k± for these T waves. 

If B ||i while both h and the chain of particles are paral¬ 
lel to z, then the dispersion relation for the coupled y and 
z waves are obtained as solutions to the coupled equations 


Poy 

POz 


-2a^ 

3d3 


2 ^yyPoy 
L^_l 


2 

- tyzPOz 

e_L 


cos(fed). 


-2a3 

3d3 



2 

- tzzPOz 

ej_ 


cos(fcd). 


( 6 ) 


These coupled y and z modes are uncoupled from the x 
mode. 

If we solve this pair of equations for poy and poz for a 
given value of k, we obtain a nonzero solution only if the 
determinant of the matrix of coefficients vanishes. For a 
given real frequency w, there will, in general, be two so¬ 
lutions for k{ui) which decay in the +z direction. These 
correspond to two branches of propagating plasmon (or 
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Fig. 2: Real and imaginary parts of 9d, the rotation angle per 
interparticle spacing (in radians), as a function of frequency, 
assuming ujcli^p = 3.5 x 10“®. Nearly horizontal chains of blue 
(+) and green (x) symbols: Drude function for metal particles 
with no damping (r —>■ oo). Black triangles and red circles: 
finite damping (ujpT = 100). In both cases, the magnetic field 
and the director of the NLC are assumed parallel to the chain 
axis, as in Fig. [T] The dotted lines merely connect the points. 


plasmon polariton) waves, with dispersion relations which 
we may write as k±[(jj). The frequency dependence ap¬ 
pears in tyz, tyy, aud tzz, which depend on w [through 
e{uj) and A(w)]. However, unlike the case where the mag¬ 
netic field is parallel to the z axis, the waves are elliptically 
rather than circularly polarized. 

Cholesteric Liquid Crystal. We now consider immers¬ 
ing the chain of metallic nanoparticles in a CLC in the 
presence of a static magnetic field with Bjjz and the chain. 
A CLC can be thought of as an NLC whose director axis 
lies perpendicular to a rotation axis (which we take to be 
z ), and which spirals about that axis with a pitch angle a 
per interparticle spacing. For a CLC, if we include only in¬ 
teractions between nearest-neighbor dipoles, the coupled 
dipole equation [Eq. ([3])] takes the form 

Pn = -^^[R.“^(zi)tC? -p^+i + R(zi)t5 -Pn-l], (7) 

as is shown in Refs. [S] and El Here = R„(z)p„ and 
R„(z) is a 2 X 2 rotation matrix for the director h{z). 
If the magnetic field lies along the z axis then the two 
T branches are coupled. One can write a 2 x 2 matrix 
equation for the coupled dipole equations in the rotated x 
and y directions. This equation is found to be 

2^3 ., 

Po =-^M(fc,u;) • po, (8) 

where po is the rotated two-component column vector 
whose components are p^o and Pyo- The components of 


the matrix M(fc,a;) are found to be 

Mxx = Cl [txx cos{kd) cos(ad) + itxy sui{kd) sin(Q!d)] 
Myy = e2[tyy cos{kd) cos{ad) + itxy sui{kd) sm{ad)] 
Mxy = e2[txy cos(kd) cos{ad) — ityy sui{kd) sm{ad)] 
Myx = ei[zta;x sin(fc(i) sin(ad) — txy cos(fc(i) cos(a(i)]. 

(9) 

1/2 

where ei = and e? = , ^ One can now determine 

the dispersion relation for the two T waves as non-trivial 
solutions to the secular equation formed from Eqs. m 
and (El). 

The most interesting result emerging from Eqs. m 
and (E)) is that the presence of a magnetic field will 
lead to dispersion relations which are non-reciprocal^ i. e., 
u}[k) /: u}{—k) in general. The magnetic field appears only 
in the off-diagonal elements txy and tyx^ which are linear 
in the field except for very large fields. The terms involv¬ 
ing txy and tyx in Eq. (E]) are multiplied by sin (fed) and 
thus change sign when k changes sign. Thus, the secu¬ 
lar equation determining uj{k) is not even in k, implying 
that the dispersion relations are non-reciprocal. The non¬ 
reciprocal nature of the dispersion relations disappears at 
B = 0 even though the off diagonal terms of M(fc,w) 
are still nonzero, because sm(kd) appears only to second 
order. Also, when the host dielectric is an NLC, the non¬ 
reciprocity vanishes because the rotation angle a = 0 and 
all terms proportional to shi(kd) vanish, even at finite B. 

In a finite magnetic field, we define the difference in 
wave number between a right-moving or left-moving wave 
as 

Aki{u}) = |Re(A:i,L)| - |Re(A:i,fl)|, (10) 

where i = 1, 2 for the two elliptical polarizations and L, R 
for either the left-moving or right-moving branch. If, for 
example, Ak{u}) ^ 0 then the left- and right-moving waves 
have different wave numbers for a given frequency and are 
non-reciprocal. 

Faraday Rotation and Ellipticity. By solving for k{uj) 
using either Eqs. ([5]) or (E)) for an NLC, or (|S|) for a CLC, 
one finds that the two modes polarized perpendicular to 
B and propagating along the nanoparticle chain have, in 
general, different wave vectors. For the NLC, we denote 
these fc_|_(a;) and fc_(a;), for the same frequency uj, while 
for the CLC, we denote them ki{u!) and k2{uj). 

We first discuss the case of an NLC host and B||z. Then 
the two solutions represent left- and right-circularly po¬ 
larized waves propagating along the chain. A linearly 
polarized mode therefore represents an equal-amplitude 
mixture of the two circularly polarized modes. This mix¬ 
ture undergoes a rotation of the plane of polarization as it 
propagates down the chain and is analogous to the usual 
Faraday effect in a bulk dielectric. The angle of rotation 
per unit chain length may be written 

= ( 11 ) 
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Fig. 3: Dispersion relations for plasmon waves elliptically po¬ 
larized in the yz plane and propagating along a chain of metal¬ 
lic nanoparticles, assuming no damping (red open squares). 
The chain is assumed immersed in an NLC with director par¬ 
allel to the chain {z) with B|| 2 ; and cUc/ojp = 3.5 x 10“®. Blue 
x’s: same as red open squares, but assuming single-particle 
damping corresponding to uipT — 100. For cup = 1.0 x 10^® 
sec“^, the chosen Ucji-Op corresponds to about 2 Tesla. 

In the absence of damping, 6 is real. If r is finite, the 
electrons in each metal particle will experience damping 
within each particle, leading to an exponential decay of 
the plasmonic waves propagating along the chain. This 
damping is automatically included in the above formal¬ 
ism, and can be seen most easily if only nearest neighbor 
coupling is included. The quantity 

0{uj) = Oi{ijj) + i02{oj) ( 12 ) 

is then the complex angle of rotation per unit length of 
a linearly polarized wave propagating along the chain of 
metal particles. By analogy with the interpretation of a 
complex 0 in a homogeneous bulk material, Re[0(aj)] repre¬ 
sents the angle of rotation of a linearly polarized wave (per 
unit length of chain), and lm[0(w)] is the corresponding 
Faraday ellipticity i. e., the amount by which the initially 
linearly polarized wave becomes elliptically polarized as it 
propagates along the chain. 

In the case of a CLC host, neither of the two T modes is 
circularly polarized in general. Thus, the propagation of a 
linearly polarized wave along the chain cannot be simply 
interpreted in terms of Faraday rotation. 

Numerical Illustrations. — We now numerically 
calculate the dispersion relations presented in the previ¬ 
ous section using the liquid crystal known as E7. This liq¬ 
uid crystal was described in experiments by Muller |29j . 
from which we take the dielectric constants ey and e_L. 
We first consider the case of an a NLC host with both 
the director and an applied magnetic field parallel to the 
metallic chain axis z. As noted earlier, the L modes are 
unaffected by a magnetic field, but the T modes are split 
into left- and right-circularly polarized waves. To illus¬ 
trate the predictions of our simple expressions, we take 


a/d = 1/3, and we assume a magnetic field such that the 
ratio Wc/wp = 3.5 x 10“®. For a typical metallic plasma 
frequency of ~ 10^® sec“^, this ratio would correspond 
to a magnetic induction B ^ 2T. We consider both the 
undamped and damped cases; in the latter, we choose 
ujpT = 100. For propagating waves we choose solutions for 
which Im[/c±] > 0 so that these waves decay to zero, as 
expected, when z ^ oo. 

The calculated dispersion relations for the two circular 
polarizations of plasmonic wave are shown in Fig. [T]with 
and without single-particle damping. The splitting be¬ 
tween the two circularly polarized T waves is not visible 
on the scale of the figure. In this, and all subsequent plots, 
we have calculated far more points than are shown in the 
Figure, so that effectively the entire range 0 < fed < tt is 
included. 

In Fig. [21 we plot the corresponding quantity 9d{k), 
the rotation angle for a distance equal to one interparti¬ 
cle spacing, in Fig. |21 When there is no damping, we find 
that both the real and imaginary parts of 9d are extremely 
small. They become larger only when damping is included, 
as we do here by setting ojpT = 100. In this case, neither 
Re[d(u;)(i] nor Im[d(w)d] exceed about 1 x 10“® radians, 
showing that a linear incident wave acquires little elliptic¬ 
ity over such distances. Since both theory and experiment 
suggest that the wave intensity typically has an exponen¬ 
tial decay length of no more than around 20 inter-particle 
spacings in realistic metallic chains, the likely Faraday ro¬ 
tation of such a wave in practice will probably not exceed a 
millidegree or two, at most, even for a field as large as 2T. 
Thus, while the rotation found here may be measurable, it 
may not be large, at least for this simple chain geometry 
with one particle per unit cell. The present expressions 
also indicate that 9d is very nearly linear in R, so a larger 
rotation could be attained by increasing B. 

For a metallic chain in a NLC where R T i, we find, 
using the same parameters and requirements as the pre¬ 
vious case, that the two non-degenerate waves (one an L 
and the other a T wave) become mixed when B ^ 0. The 
dispersion relations, again with and without damping, are 
plotted in Fig. [31 When compared to previous work in 
Ref. ini the dispersion relations in Fig. [3] are modihed be¬ 
cause of the finite damping and presence of the magnetic 
field. However, the change produced by the magnetic field 
in these dispersion relations is not visible in the figures. 

Finally, we discuss the case of metallic chain parallel to 
the z axis, subjected to a magnetic field along the z axis, 
and immersed in a CLC whose twist axis is also parallel 
to z. Using the same host dielectric constants given above 
and a twist angle of ad = 7r/6, we show in Fig. SKa) the 
resulting dispersion relations, i. e., ui/tCp plotted against 
I fed I, for the two transverse branches. In particular, we 
show both transverse branches for a right-moving wave 
(black symbols) and left-moving wave (red symbols) giv¬ 
ing a total of 4 plots shown in Fig. SKa). The separation 
between the two T branches is on the order of 0.05 uj/uip 
for allfc. In Fig. HKb), we plot the corresponding difference 
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Fig. 4: (a). Black + symbols: the two dispersion relations for right-moving transverse plasmon waves propagating along a chain 
of metallic nanoparticles immersed in a CLC host with B||£, plotted as a function of \kd\. Red x’s: same quantities but for 
left-moving plasmon waves. We assume that = 3.5 x 10“®, uipT = 100, and ad = tv/6. (b). The difference between the 

quantities |A:id|, i. e., \Akid\ = \Re{ki^Ld)\ — |i?e(fci,_Rci)|, for left-propagating and right-propagating modes, as given by Eq. (1101) . 
for the two elliptical polarizations. L, R refer to the left-moving or right-moving waves, and i (i = 1,2) labels the two branches 
for each direction. Note that a non-zero value of Aki{uj)d implies that for a given frequency the left- and right-traveling waves 
have different wave vectors. Numerically, we find that Aki{ij)d is independent of i. 


in wavenumber between the left- and right-moving waves 
as Aki{uj)d. For the parameters used, this difference turns 
out to be the same, to within numerical accuracy, for both 
the transverse branches, and hence we show only a sin¬ 
gle plot. Since Aki(u})d is nonzero in a wide frequency 
range, the wave propagation is indeed non-reciprocal in 
this range. One-way wave propagation may occur in part 
of this range. Such propagation would occur when, at 
particular frequencies, waves can propagate only in one of 
the two directions. We have not, however, checked this 
possibility numerically. 

Discussion. — The present numerical calculations 
omit several potentially important factors which could al¬ 
ter the numerical results. Among these are the effects of 
particles beyond the nearest neighbors on the dispersion 
relations m, the (possibly large) influence of the par¬ 
ticles in disrupting the director orientation of the liquid 
crystalline host , and the effects of radiative damp¬ 

ing m on the dispersion relations. Nonetheless, we be¬ 
lieve that our calculations correctly describe, at least qual¬ 
itatively, how the surface plasmon dispersion relations are 
affected by the combined influence of a liquid crystalline 
host and an applied magnetic field. 

It should be noted that the magnetic field effects de¬ 
scribed in this paper are numerically very small, for the pa¬ 
rameters investigated. The smallness is caused mainly by 
the small value of the ratio Wc/wp, taken here as 3.5 x 10“^. 
To increase this ratio, one could either increase Wc (by 


raising the magnetic field strength), or decrease tOp (by re¬ 
ducing the free carrier density in the metal particle). We 
speculate this latter effect could be achieved by studying 
n-type semiconductor nanoparticles, if these could be pre¬ 
pared in a free carrier regime without having to consider 
quantized energy levels. The same calculations should also 
be carried out for more realistic metallic dielectric func¬ 
tions than the Drude model. 

In summary, we have calculated the dispersion relations 
for plasmonic waves propagating along a chain of metallic 
nanoparticles immersed in a liquid crystal and subjected 
to an applied magnetic field. For a magnetic field paral¬ 
lel to the chain and director axis of the NLC, a linearly 
polarized wave is Faraday-rotated by an amount propor¬ 
tional to the magnetic field strength. For a CLC host and 
a magnetic field parallel to the chain, the transverse wave 
solutions become non-reciprocal (left and right-traveling 
waves having different dispersion relations) and there may 
be frequency ranges in which waves can propagate only in 
one direction. Thus, plasmonic wave propagation can be 
tuned, either by a liquid crystalline host or a magnetic 
field, or both. In the future, it may be possible to detect 
some of these effects in experiments, and to use some of 
the predicted properties for applications, e. g., in optical 
circuit design. 
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